Nonequilibrium steady states in contact: 
Approximate thermodynamic structure and zero-th law for driven lattice gases 
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We explore driven lattice gases for the existence of an intensive thermodynamic variable which could 
determine "equilibration" between two nonequilibrium steady-state systems kept in weak contact. 
In simulations, we find that these systems satisfy surprisingly simple thermodynamic laws, such as 
the zero-th law and the fluctuation-response relation between the particle-number fluctuation and 
the corresponding susceptibility remarkably well. However at higher densities, small but observable 
deviations from these laws occur due to nontrivial contact dynamics and the presence of long-range 
spatial correlations. 
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Among the wide class of nonequilibrium systems, an 
important and ubiquitous subclass are those which have 
a nonequilibrium steady state (NESS). Unlike in equilib- 
rium, a system in a NESS has a steady current but its 
macroscopic properties, like in equilibrium, are still in- 
dependent of time. In contrast to equilibrium systems, 
there is no well founded thermodynamic theory even for 
this conceptually simplest class of nonequilibrium sys- 
tems. Intensive studies attempting to construct a suit- 
able statistical mechanical framework where macroscopic 
properties and thermodynamic states may be character- 
ized in a simple way have not yet converged to a universal 
picture dHU. 

At the heart of equilibrium thermodynamics is the 
zero-th law which is a consequence of equalization of in- 
tensive thermodynamic variables when two systems are 
in contact. For example, when two systems with the same 
temperature are allowed to exchange particles with the 
total number of particles conserved, the final equilibrium 
state is determined by equalization of the chemical po- 
tentials of the two, obtained by minimizing the total free 
energy. For NESSs, we ask the same: What happens if 
two NESSs are brought into contact? 

Recently, there have been attempts to define an inten- 
sive thermodynamic variable for systems such as driven 
granular systems 0] , static granular assemblies of blocked 
states formed by weak driving Q and a class of exactly 
solvable models motivated by inelastic granular collisions 
Q. More generally, there has been a prescription to de- 
fine such a variable for systems in NESSs by invoking a 
hypothesis, called the asymptotic factorization property, 
which has been shown to be satisfied for a class of systems 
having short-range spatial correlations Q . 

For driven diffusive systems like the paradigmatic 
stochastic lattice gases [h], E3, which have long-range 
spatial correlations, the situation is less clear. Previ- 
ously, motivated by equilibrium thermodynamics which 
has a rigorous basis in terms of the large-deviation princi- 
ple (LDP) [l3j . a hypothesis of the existence of LDP has 
been put forward for these systems but not yet rig- 



orously established. By operationally defining a pressure 
and a chemical potential, a numerical study |l4j indi- 
cates that a Maxwell relation is satisfied and there may 
indeed exist a large-deviation function analogous to the 
equilibrium free energy. 

In equilibrium, the existence of an intensive variable 
hinges crucially on the local thermodynamic properties 
of a system, i.e., if the system is divided into subsys- 
tems large compared to the microscopic scales, the fluc- 
tuations in the individual subsystems are independent of 
each other as a consequence of the short-range spatial 
correlations in the system. In contrast, the driven sys- 
tems have generic long-range spatial correlations fl5| . In 
this situation, it is not obvious that the system could be 
divided into independent subsystems and intensive vari- 
ables analogous to those in equilibrium could be defined. 

In this paper we explore by simulations the "equili- 
bration" between two driven lattice gases upon contact. 
Interestingly, we find that, to a very good approximation, 
there is an intensive variable, like equilibrium chemical 
potential, which determines the final steady state while 
two such driven systems are allowed to exchange parti- 
cles. Concomitantly, the zero-th law of thermodynamics 
is satisfied remarkably well. Moreover, a fluctuation- 
response relation between the fluctuations in particle- 
number and the corresponding susceptibility is also well 
satisfied. However, at higher densities, there are small 
but observable deviations from these simple thermody- 
namic laws due to nontrivial contact dynamics and the 
presence of long-range spatial correlations. 

We consider two systems of volume V\ and Va, con- 
nected at a finite set of points V\ and V 2 which are sub- 
sets of Vi and V 2 respectively, with Vi , V 2 -C V\ , V2 (see 
Fig. [IJ. The two systems can interact and exchange 
particles with each other only at the contact. The en- 
ergy H of the two systems combined is given by H = 

KiEv(riHn') + K2Ev(r 2 )v(r 2 / ) + #E*7(ri)??(r 2 ) 
where sums are over nearest-neighbor sites with ri, iV 6 
Vi, ?i € Vi and r 2 ,r 2 ' £ V 2 , ? 2 & V 2 . A site r can 
be occupied by at most one particle and the occupation 
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FIG. 1: A schematic diagram of two nonequilibrium steady 
states with contact region Vi and V2. 



variable n(r) is 1 or if the site is occupied or unoc- 
cupied, respectively. K\, K2 and K are the interaction 
strengths among particles for system 1 and 2, and at 
the contact, respectively. In the simulations, we con- 
sider two-dimensional systems (V = L x L) with peri- 
odic boundaries in both directions. We choose the jump 
rate w(C'\C) from a configuration C to C according 



to the local detailed balance condition [121] : the jump 
rate from a site r to its unoccupied nearest neighbor r' 
obeys w{C'\C) = w(C\C) cxp[-AH + E(x' - a:)] where 
AH = H(C') - H(C), E is the driving field along the 
x-direction, and x and x' are x-components of r and r' 
(ksT = 1, ks the Boltzmann constant, T temperature). 

We choose E = E\ when r, r' G V\, E = E 2 when 
r, r' G V 2 and E = otherwise. There is no driving 
field along the bonds connecting the two systems. We 
choose K\ , K2 > and E\ ,E> so that systems are in 
the disordered (fluid) phase [lOj . For E\ = E2 = 0, the 
combined system has the equilibrium Boltzmann distri- 
bution ~ exp[— H(C)]. For Ei,E 2 / 0, there are currents 
in the steady states and the steady state distribution is, 
in general, unknown. Driven bilayer systems were stud- 
ied previously [U| where particles jump from one layer 
to the other at any site as opposed to the case here with 
possibility of particle transfer only at a small contact 
area. 

We first report an a priori surprising observation sug- 
gesting an effective zero-th law for systems in NESS. 
When two systems are brought into contact, after relax- 
ation involving exchange of particles, in the final "equi- 
librated" steady state there is no net current across the 
contact region. This property allows an operational def- 
inition of a chemical potential of a NESS as follows. A 
driven system is brought into contact with an equilib- 
rium system EQ whose chemical potential fi is known as 
a function of density. In the final steady state, the chem- 
ical potential fi of EQ is assigned to the driven system. 
We choose a system of noninteracting hard-core parti- 
cles as EQ with density being no and /i = —(ds/dno) = 
ln[no/(l — no)] with s = — [no In no + (1 — no) ln(l — no)] 
the equilibrium entropy per lattice site. By varying no 
of EQ in contact with a NESS in consideration, one can 
get the density versus chemical potential curve for the 
NESS as shown in Fig. [5] (bottom panel). The surpris- 



ing observation is that, if two NESSs, NESSi and NESS 2 
(chosen such that they have operationally the same n 
but different densities), are brought together, the respec- 
tive densities do not change upon contact. Moreover, if 
we bring together two NESSs with the same density but 
different /it, particles will flow from the higher to lower 
chemical potential till the respective densities correspond 
to the same [i as indicated with arrows in Fig. [2] (bottom 
panel) . 

Thus, if two systems are separately equilibrated with 
a common system with a fixed density, they will also be 
equilibrated amongst themselves. Consider, e.g., two sys- 
tems NESSi and NESS2 kept in contact and having two 
equilibrated density profiles, with density ni and density 
n 2l respectively. Then, a third system EQi is separately 
brought into contact with NESS2 and the density of EQi 
is tuned to such that NESS2 keeps its density n 2 un- 
changed in the equilibrated state. Now, if NESSi with 
density n\ and EQi with density n^, are brought into con- 
tact, the two density profiles remain almost unchanged, 
confirming the zero-th law (see Fig. [2j top panel and the 
explanations in the caption). 

These systems are indeed far away from equilibrium 
since the numerical values of the currents in NESSi and 
NESS2 in the bottom panel of Fig. [2] are approximately 
2/3 and 1/3 of the respective maximum currents. Like- 
wise, in the top panel, NESS2 with density ss 0.5 has a 
homogeneous disordered state in contrast to the corre- 
sponding equilibrium system, with K = 2, E = and 
the same density, which has a symmetry-broken phase 
with different sub-lattice densities (lpj . 

The chemical potential of a system may also be mea- 
sured by keeping it in contact with any other equilibrium 
system, not necessarily a noninteracting hardcore one. 
In Fig. [3] we have plotted densities versus chemical po- 
tentials for a system separately in contact with systems 
with different contact area as well as nonzero interaction 
strength K 1 the result being in good agreement with the 
zero-th law. 

The existence of a zero-th law would be a consequence 
of a putative large-deviation principle (LDP) [13j. To 
elucidate it briefly, let us consider two systems which 
can exchange particles such that N\ + N 2 = N = const , 
N\, N2 the number of particles in systems 1 and 2 re- 
spectively. Assuming that the LDP holds, the probabil- 
ity P(Ni,N 2 ) of a large deviation in N±,N 2 is given by 



P(N!,N 2 ) 



a VlSi(rti) _,V2S2("2) 



] e" 5(Ar) in the limit of 



Ni,N 2 ,Vi,V2 > 1 where ni = N1/V1 and n 2 = N 2 /V 2 
being finite and exp[— S'(iV)] the normalization constant 
( £ ~' implying equality in terms of logarithm). The func- 
tions S\(ni), S2(n2) are called the large-deviation func- 
tions (LDF). In writing so, the correlation between sys- 
tems has been neglected as a boundary-effect in the limit 
of large volume. This assumption of a product measure of 
P(Ni,N 2 ) essentially implies that the LDFs si(ni) and 
s 2 (n 2 ) are local function of the respective densities. The 
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FIG. 2: Bottom panel: The plot of density n vs. chemical 
potential fi for two 120 x 120 systems, NESSi with K = 1, 
£ = 2 (circles) and NESS 2 with A = E = 2 (triangles) with 
A' = 0. Arrows indicate how the density changes if NESSi 
and NESS2 with the same initial density n = 0.30 (denoted 
by middle horizontal line) are brought into contact, reaching 
respective final densities n ~ 0.33 (denoted by top horizontal 
line) and n ~ 0.27 (denoted by bottom horizontal line) with 
equal chemical potential /1 ~ 0.36. Top panel: Numerical 
experiments to test zero-th law (K = 0) - (1) NESSi (A = 
4, E — 6, L = 120) with density n\ (bottom red profile) 
equilibrated with NESS2 (A = E = 2, L = 110) with density 
712 (middle blue profile), (2) NESS2 with density 712 (middle 
magenta profile) equilibrated with EQi (A = 1, E = 0, L = 
100) with density 713 (top grey profile), and (3) NESSi with 
density n'i (bottom green profile) equilibrated with EQi with 
density n' 3 (top black profile) where n! x « ri\ and n' 3 ~ 113. 




FIG. 3: Density n vs. chemical potential fi is plotted for a 
120 x 120 NESSi (A = 1, E = 2) for the following cases. 
Bottom panel: NESSi separately in contact with 120 x 120 
systems with EQi (A = E = 0, red) and EQ 2 (A = 1, E = 0, 
blue). Middle panel: NESSi separately in contact with EQi 
(blue) and EQ2 (black) with 4x1 contact area. Top panel: 
NESSi separately in contact with EQi (blue) and EQ2 (black) 
for nonzero interaction strength K = 1. In middle and top 
panel, n vs. fi plot is compared with that obtained for NESSi 
in contact with EQi with 2x2 contact area. 



macroscopic state, under the constraint N1+N2 = const, 
is determined by maximizing In P (iVi,^) where the 
chemical potentials fii = —dsi/dni and u 2 = —ds2/dri2 
being equal in the final steady state. Clearly the conse- 
quence of the LDP is a zero-th law as presented in Fig. 
[21 Another interesting consequence of a putative LDP 
would be a relation between the susceptibility and the 
fluctuation in particle-number of the system 1 in a NESS 
when it is in contact with the system 2 being a large 
reservoir characterized by a chemical potential u. Then, 
one gets the following fluctuation-response relation as in 
equilibrium, 



X = ^ = ((Nf) - (iV,) 2 ) - 



'iVi 



(1) 



We first proceed to test this relation for a NESS 
in contact with an equilibrium reservoir with density 
no, consisting of noninteracting hardcore particles with 
u = —(ds/dno) = ln[no/(l — "o)]- For better nu- 
merical accuracy, we check the integrated version of 
Eq. d]by defining the integrated susceptibility I X ((J.) = 
^{d{N x )/d(t)dn = (NiQi)) - (Ntifio)) and the into- 
grated fluctuation I a (fJ.) = J ^(cr^^dfj,. 

We take a 20 x 20 nonequilibrium system NESSi with 
K = 1, E = 2 and keep it in contact with a 100 x 100 
equilibrium reservoir of noninteracting hardcore particles 
(K = E = 0), called RESi. Then we vary the chemical 
potential u (or equivalently the density no) of RESi in 
small steps from an initial value uq = —3.5 and calculate 
a% i for each value of n. We repeat this procedure by 
keeping NESSi separately in contact with various other 
reservoirs of size 100 x 100 whose chemical potentials can 
be measured by keeping these reservoirs in contact with 
the RESi. In Fig. |H we plot l x {p) and I a {p) as a func- 
tion of [i. Provided that the Eq. [T]is valid, all the curves 
should fall on each other. Up to chemical potential u ~ 1, 
we observe a quite good collapse within the numerical ac- 
curacy. We also consider two different systems, NESSi 
with K = 1, E = 2 and NESS 2 with K = E = 2, sepa- 
rately in contact with RESi. In the inset of Fig. 21 we 
plot I x (n) and I a {n) which are in good agreement with 
the fluctuation relation in Eq. [TJ 

At higher chemical potentials, there are observable de- 
viations from this simple thermodynamic behavior. In 
Fig. HI the I x vs. u and I a vs. fj, curves do not 
fall on each other for fj, > 1. Correspondingly, in this 
density regime, the zero-th law does not hold strictly 
as seen in Fig. |3l However, these violations are not 
simply due to a finite-size effect and persist for much 
larger system sizes. To investigate the possible rea- 
sons for the violations, we also study the behavior of 
spatial density correlation functions for various densi- 
ties. Unlike equilibrium systems, the nonequilibrium 
systems, due to the presence of a driving field, are ex- 
pected to have generic long-ranged spatial correlations 
\m 15-171 , arising because the structure factor S(q x -> Qy); 
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FIG. 4: Integrated susceptibilities I x (squares) and fluctua- 
tions I a (circles) vs. chemical potential fj, are plotted for a 
20 x 20 NESSi with K = 1, E = 2, separately in contact 
with five different 100 x 100 reservoirs: (1) K = E = (red), 
(2) K — 1. E — 2 (green), (3) = 1, E = 4 (blue), (4) 
K = 2, E = 4 (sky-blue), (5) AT = 2, S = 6 (magenta). Inset: 
Same quantities are plotted for two different 20 x 20 systems, 
NESSi with K = 1, E = 2 (blue) and NESS 2 K = E = 2 
(red), in contact with a reservoir (A' = E — 0, L = 100). 



i.e., the Fourier transform of the spatial density correla- 
tion function, becomes singular when q Xl q y —> with 
R = [lim^o 5(0, q y )}/[\im q ^ S(q x , 0)] ^ 1. This gives 
rise to the long-range spatial correlations decaying as 
A/r d with distance r in d dimension with the amplitude 
A tx (R - 1) [lot. In Fig. we have plotted the ratio 
i? versus density n and the structure factors S(q x , 0) and 
5(0, q y ) for two different densities in inset of Fig. [5] for a 
NESS with K = 1, E = 2, L = 120. The ratio R deviates 
from 1 more strongly with increasing density, indicat- 
ing an increase of long-range correlations which manifest 
themselves through the nonlocal effect of the contact. In 
Figs. [3]and|H collapse of various n vs. /i and I a vs. \i (or 
I x vs. fi) curves for different contact dynamics are not 
very good for \i > 1. Clearly, the effect of the contact is 
felt throughout the systems and this nontrivially changes 
the corresponding thermodynamic properties. Therefore, 
the break-down of an exact equilibrium-like structure or, 
in other words the break-down of the product-measure 
assumption in the LDP, indicates the important role of 
the contact dynamics [9| and the long-range correlations 
in a driven system. 

In summary, our numerical study of coupled driven 
lattice gases has revealed a surprisingly simple thermo- 
dynamic structure with an effective zero-th law like be- 
havior concerning exchange of particles and the corre- 
sponding fluctuation-response relation. This thermody- 
namic structure is not exact since there are small but 
observable deviations at higher densities. Their physi- 
cal origin is rooted in the nontrivial contact dynamics 
and the presence of long-range spatial correlations which 
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FIG. 5: The ratio R = [lim^o S{0, q y )]/[lim qx ^ S(q x , 0)] 
vs. density n is plotted for a NESS with K = 1, E — 2, L — 
120. Inset: The structure factors S(q x ,0) (open points) and 
S(0, q y ) (filled points) vs. q x and q y respectively are plotted 
for densities n — 0.1 (red) and n — 0.5 (blue). 



invalidate the asymptotic factorization property. As an 
open question, it would be interesting to sec whether sys- 
tems with more than one conserved quantity (e.g., mod- 
els with two species) exhibit a similar behavior. Finally 
our study prompts the question whether such an approx- 
imate thermodynamic structure is typical just for driven 
lattice gases or generically occurs in other coupled NESSs 
with long-range correlations as well. 

We thank J. Krug and R. K. P. Zia for discussions. 
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